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Abstract 

Recently, the first-named author gave a classification of 3D consistent 6-tuples 
of quad-equations with the tetrahedron property; several novel asymmetric 6-tuplcs 
have been found. Due to 3D consistency these 6-tuples can be extended to discrete 
integrable systems on Z m . We establish Lagrangian structures and flip-invariance of 
the action functional for the class of discrete integrable systems involving equations 
for which some of the biquadratics are non-degenerate and some are degenerate. 
This class covers, among others, some of the above mentioned novel systems. 



1. Introduction 



One of the definitions of integrability of discrete equations, which becomes increasingly 
popular in the recent years, is based on the notion of multidimensional consistency. 
For two-dimensional systems, this notion was clearly formulated first in [NW01J, and 
it was pushed forward as a synonym of integrability in |BS02| Nij02| . The outstanding 
importance of 3D consistency in the theory of discrete integrable systems became evident 
no later than with the appearance of the well-known ABS-classification of integrable 
equations in [ABS03J. This classification deals with quad-equations which are set on 
the faces of an elementary cube in such a manner that all faces carry similar equations 
differing only by the parameter values assigned to the edges of a cube. Moreover, each 
single equation admits a symmetry group. Thus, ABS-equations from [ABS03J can 
be extended in a straightforward manner to the whole of Z m . In |ABS09| . a relaxed 
definition of 3D consistency was introduced: the faces of an elementary cube are allowed 
to carry a priori different quad-equations. The classification performed in |IABS09| is 
restricted to so-called equations of type Q, i.e., those equations whose biquadratics are 
all non-degenerate (a precise definition will be recalled in Section [2]). In [AtkO§J and 
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2. Quad-equations on a single square 



[BolSlOj, numerous asymmetric systems of quad-equations have been studied, which also 
include equations of type H (with degenerate biquadratics). A classification of such 
systems has been given in |Bolllj . This classification covers the majority of systems 
from [Atk08j and all systems from jBolSlOj . as well as equations from |Hie04] . |LY09| 
and |HV10] , Moreover, it contains also a number of novel systems. The results of this 
local classification lead to integrable lattice systems via a procedure of reflecting the 
cubes in a suitable way (sec [BolllJ for details). As shown in [Atk08j, [BolSlOj, |Bolll| . 
the asymmetric systems still can be seen as families of Backlund transformations, and 
they lead to zero curvature representations of participating quad-equations. 

After the pioneering work |MV91j the variational (Lagrangian) structure of discrete 
integrable systems is a topic which receives more and more attention. A Lagrangian 
formulation of systems from the ABS list is given in [ABS03J. The new idea in [LN09J 
was to extend the action functional of |ABS03| to a multidimensional lattice, which is 
meaningful due to the invariance of the action under elementary 3D flips of 2D quad- 
surfaces in Z m . This remarkable and fundamental flip invariance property was proven in 
|LN09] for some particular cases from the ABS list and in |BS10| for the whole list. 

In the present paper, we introduce a Lagrangian formulation for an important subclass 
of systems from the classification in |Bolllj , namely those which include type H equations 
not all of whose biquadratics are degenerate, and prove the flip invariance of action 
functionals for these systems. 

The outline of the paper is as follows: In Sections [2] and [3] we introduce the necessary 
notions and notation and recapitulate the classification results from |Bolll| which are 
relevant for the present paper. Then, in Section [4] we use the so-called three-leg form 
of quad-equations to construct the Lagrangians for our equations. The fundamental 
property of Lagrangians on an elementary quadrilateral constitutes one of the main 



results of the paper (Theorem 4.4). The Lagrangians are elementary building blocks of 
the action functional constructed in Section [5] which also contains the second main result 
of the paper, namely the flip invariance of the action functional for systems of our class. 

2. Quad-equations on a single square 

We start with introducing relevant objects and notations. A quad- equation is a relation 
of the type Q (xi, x 2 , £3, £4) = 0, where Q £ C [x±, x 2 , £3, X4] is an irreducible multi-affine 
polynomial. It is convenient to visualize a quad-equation by an elementary square whose 
vertices carry the four fields x±, x 2 , %3, X4, cf. Figure [T] For quad-equations without 
pre-supposed symmetries, the natural classification problem is posed modulo the action 
of four independent Mobius transformations of all four variables: 
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In the simplest situation, a quad-equation is thought of as an elementary building block 
of a discrete system Q{x m ^ n , x m +i in , 2Cm+l,n+l> x m ,n+i) = for a function x : 1? — > C, 
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2. Quad-equations on a single square 



but we will also encounter more tricky ways of composing discrete systems on Z 2 from 
different quad-equations. 

A very useful tool for studying and classifying quad-equations, or, better, multi-affine 
polynomials of four variables, are their biquadratics. For any two (different) indices 
i,j 6 {1, 2, 3, 4}, we define the corresponding biquadratic as follows: 



Q ' — Q (xi, Xj) — QxkQxi QQxk 



X( 1 



where {A;, £} is the complement of {i, j} in {1, 2, 3, 4}. A biquadratic polynomial h (x, y) 
is called non- degenerate if no polynomial in its equivalence class with respect to Mobius 
transformations in x and y is divisible by a factor x — a\ or y — ct<i (with ttj G C). Oth- 
erwise, h is called degenerate. Here, of course, the action of two Mobius transformations 
on h is defined as 

t / s i , , s2/ . , x2 u (aix + bi a 2 y + b 2 \ 

h(x,y) ~» (cix + dx) {c 2 y + d 2 ) ■ h[ — , — — J . 

Vcix + di c 2 y + d 2 J 

We recall also that an important role in studying the multi-affine polynomials Q belongs 
to the four discriminants 

,•'(,•;) «l;:{r 2Q';i,,Q- J 

(it can be shown that the result does not depend on j). These are quartic polynomials 
which can be assigned to the vertices of an elementary square. Note that the action of a 
Mobius transformation on a quartic polynomial r(x) is defined as 

ax + 6 N 



r(x) [cx + d) -r 

\cx + d) 



cx + d' 

A multi-affine polynomial Q is said to be of type Q if all its biquadratics are non- 
degenerate and of type H otherwise. It was proven in [BolllJ that in the latter case 
either four out of six biquadratics are degenerate, and then Q is said to be of type H^, 
or all six biquadratics are degenerate, and then Q is said to be of type H 6 . 

A complete classification of multi-affine polynomials Q up to Mobius transformations 
in x\, ... ,X4 is given in [BolllJ. In the present paper, we will be only interested in 3D 
consistent six-tuples of quad-equations containing at least one equation of type H 4 but no 
equations of type H 6 . (Recall that six-tuples consisting solely of type Q equations were 
already classified in [ABS09J, and the flip invariance of the action functional for such 
systems was completely discussed in [BSlOj.) As shown in |Bolll| . the list of canonical 
normal forms of polynomials of type H 4 consists of three items (which actually appeared 
already in |ABS09j ): 

H{: Q = (xi- x 3 ) (x 2 - X4) + (a 2 - a{) (l - e 2 x 2 x 4 ), 

H| : Q = i x i - ^3) {x 2 - X4) + (a 2 - a{) (a?i + x 2 + x 3 + X4) + a 2 - a\ 

£ go 

-- {a 2 - a\) (2x 2 + cti + a 2 ) (2rr 4 + a\ + a 2 ) — - (a 2 - a>i) , 



m-. Q = e 2ai (xix 2 + X3X4) - e 2a2 (xtx 4 + x 2 x 3 ) + (e 4ai - e 4a2 ) ( 5 2 + 



t 2 X2X4 

+012) 
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2. Quad-equations on a single square 



Thus, the members of this list, like those of the list Q, depend on two parameters ai, a 2 
which can be assigned to the pairs of opposite edges of the elementary square. However, 
unlike the Q case, where all vertices are on the equal footing (which can be visualized as on 
Figure la), for any type H 4 polynomial, the four discriminants assigned to the vertices 
belong to two different types (at least for e / 0). Up to constant factors K 2 (ai,a 2 ) 



depending on the parameters only and listed in Section A. 2 (cf. also Lemma 4.1), one 
has: 



H{: r 1 (xi)=e 2 , r 2 (x 2 ) = 0, r 3 (x 3 ) 



r\x 4 ) = ; 



HI: 



r 1 (xi) = 1 + 4exi, r 2 (x 2 ) = 1, ?~ 3 (x3) = 1 + 4ex3, r 4 (x4) = 1; 
r l (xi) = x 2 — 45 2 e 2 , r 2 (x 2 ) = x 2 , r 3 (x 3 ) = x 2 — 45 2 e 2 , r 4 (x4) = x\. 



This leads to the coloring of the vertices in two black and two white ones. We will think 
of the black vertices as carrying the parameter e, while for the white vertices e = 0. 

It will be useful to visually keep track of the distribution of the four degenerate and 
two non-degenerate biquadratics of a type H polynomial. We will indicate the two non- 
degenerate ones by thick lines. When composing discrete systems from quad-equations, it 
becomes crucial how to arrange the variables (fields) x\, X2, X3 and X4 over the vertices of 
an elementary square. There are two different (up to rotations of an elementary square) 
quad-equations of type H 4 , corresponding to two different biquadratics patterns, as shown 
on Figure [T] 




X\ OL\ X2 

(a) Quad-equation of type Q 





(b) Rhombic version of a 
quad-equation of type H 4 



(c) Trapezoidal version of a 
quad-equation of type H 4 



Figure 1: Biquadratics patterns of quad-equations; non-degenerate biquadratics are in- 
dicated by thick lines 



The first arrangement corresponds to assigning the two black vertices to one diagonal 
of the elementary square and the two white vertices to another diagonal, see Figure 
lb We will call these quad-equations the rhombic version of due to the symmetry 
properties 



Q (xi,x 2 , x 3 , x 4 ; ai, £22, e) = -Q (x 3 ,x 2 , X\,x& a 2 , ai, e) 
of the corresponding polynomials. 



-Q(xi,X4,x 3 ,x 2 ;a2,ai,e) 
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3. 3D consistent six-tuples of quad-equations 



The second version of a quad-equation is obtained from the canonical normal form Hi 
by the switching X2 O £3 and a\ •<-)■ «2 — As a consequence, the two black vertices 
share an edge rather than a diagonal, and the same holds true for the two white vertices. 
This version of a quad-equation possesses only one axis symmetry which interchanges the 



vertices of the same color (see Figure lc). We will call the corresponding quad-equations 
the trapezoidal version of Hi. 

By the way, the biquadratics pattern provides us with an unambiguous way of bi- 
coloring the vertices of the square even in the case e = 0, when all four discriminants 
are the same (this corresponds to the list H from [ABS03J): each of the two thick lines 
corresponding to the two non-degenerate biquadratics connects two vertices of the same 
color. 

Concerning discrete systems composed of quad-equations, we mention that the rhombic 
version of any equation Hi can be imposed on any bipartite quad-graph (cf. also [XP09J). 
As for the trapezoidal version of H^ it can be imposed on the regular square lattice 
Z 2 with horizontal (or vertical) rows colored alternately black and white, see Figure 
[2] Integrability of such a non-autonomous system of quad-equations understood as 3D 
consistency has been discussed in |Bolll| . 



a 2 



cti a 3 a 5 

Figure 2: Integrable system on Z 2 consisting of quad-equations Hi 



3. 3D consistent six-tuples of quad-equations 

We will now consider six-tuples of (a priori different) quad-equations assigned to the 
faces of a 3D cube: 

A (x, x x , x 2 , Z12) = 0, A (23, X13, X23, £123) = 0, 

B (X,X 2 ,X 3 ,X 2 3) = 0, B (X1,X12, Xi 3 ,Xi23) = 0, (1) 

C (X, Xx, X 3 , X 1S ) =0, C (X 2 ,X 12 ,X23, x 123 ) = 0, 



see Figure 3a Such a six-tuple is 3D consistent if, for arbitrary initial data x, x\, 22 
and X3, the three values for X123 (calculated by using A = 0, B = or C = 0) coincide. 
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3. 3D consistent six-tuples of quad-equations 



A 3D consistent six-tuple is said to possess the tetrahedron property if there exist two 
polynomials K and K such that the equations 

K (x, x 12 , X13, X23) = 0, K (x 1 ,x 2 , x 3 , X123) = 

are satisfied for every solution of the six-tuple. It can be shown that the polynomials K 
and K are multi-affine and irreducible (see (BplllJ). 



x 3 



B- 



X23 



A 




C 



•i'2 



.4 



B 



X123 




X12 




X123 



X12 



X13 



(a) A 3D consistent six-tuple (b) Making tetrahedra to faces 

Figure 3: Equations on a cube 

A complete classification of 3D consistent six-tuples ([I]) possessing the tetrahedron 
property, modulo (possibly different) Mobius transformations of the eight fields x, Xi, 
Xij, X123, is given in jBolllj . We mention one of the important general properties of such 
six-tuples: 

Lemma 3.1. Equations on opposite faces are of the same type and have the same bi- 
quadratics pattern. 

In the present article, we consider only six-tuples ([!]) whose equations are not all of 



type Q and which do not contain type H equations. Taking into account Lemma 3.1 



one easily sees that three different combinatorial arrangements of biquadratics patterns 
are possible, as indicated on Figure [4] Arrangements on Figures [4a] and [4b] are related 
to each other by the following construction, sec [BolllJ for a proof: 

Lemma 3.2. Consider a 3D consistent six- tuple possessing the tetrahedron property 
expressed by the two equations 



K (x,xi2,a;i3,X23) 

Then the six-tuple 

K(x, Xi2,Xi3,X23) = 
B (x,X 2 ,X3,X23) = 
C (x,Xl,X3,Xi 3 ) = 



K (x 1 ,X 2 ,X 3 ,Xi23) 



0. 



0, 

0, 
0, 



K (x 1 ,X 2 ,X 3 ,X 12 3) = 0, 
S(xi,Xl2,2;i3,Xi23) = 0, 
C(x 2 ,Xl2,a;23,^123) = 0, 



(2) 
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3. 3D consistent six-tuples of quad-equations 




(a) First case: all face 
equations of type H 4 , 
tetrahedron equations 
of type Q 




(b) Second case: two pairs 
of face equations and 
tetrahedron equations 
of type H 4 , one pair of 
face equations of type 

Q 




(c) Third case: all face 
and tetrahedron equa- 
tions of type H 4 



Figure 4: Biquadratics patterns 



assigned to the faces of a 3D cube as on Figure 3b is 3D consistent and possesses the 
tetrahedron property expressed by the two equations 

A (x, x 1 ,x 2 ,xi 2 ) = 0, A (x 3 , Xi 3 , X 2 3, x 123 ) = 0. 

3D consistency of ^ is understood as the property of the initial value problem with the 
initial date x, x%, X\% and x 23 . 

Now we can formulate the main result of jBolllj concerning 3D consistent six-tuples 
containing equations of type H 4 but no equations of type H 6 . In this theorem, we use 
normal forms Q k of equations of the list Q which are slightly different from those adopted 
in |ABS03llABltl9] . namely 

Q\ : Q = a\ (x\X2 + x 3 x 4 ) - a 2 (xix 4 + x 2 x 3 ) - {a\ - a 2 ) {x\x 3 + x 2 x 4 ) 
+e 2 a\a 2 {a\ — a 2 ), 

Q\: Q = a.\ (x\x 2 + x 3 X4)-a 2 (21X4 + x 2 x 3 )-(a>i - a 2 ) (xi 23 + x 2 Xi)+aia 2 (a\ - ct 2 ) 
+ea\a 2 {a\ - a 2 ) (xi + x 2 + x 3 + x 4 ) - e 2 a\a 2 (a\ - a 2 ) (a\ - a\a 2 + a|), 

Q|: Q = sinh (2ai) {x\x 2 + X3X4) — sinh (2a 2 ) (xix^ + x 2 x 3 ) 

— sinh (2 (ai — a 2 )) {x\x 3 + x 2 x^)— 4<5 2 e 2 sinh (2«i) sinh (2a 2 ) sinh (2 {a\ — a 2 )). 

Theorem 3.3. Modulo independent Mobius transformations of all fields, there are nine 
3D consistent six-tuples of quad- equations with the tetrahedron property which contain 
equations of type Br but no equations of type H 6 : 



a) Three systems with the biquadratics pattern as on Figure \4a\ with all six face equa- 
tions being rhombic H^, and the tetrahedron equations Qi (for k = 1,2,3). 

b) Three systems with the biquadratics pattern as on Figure with two pairs of face 
equations being trapezoidal H^, one pair of face equations Q e k , and the tetrahedron 
equations H k (for k = 1,2,3). 
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4. Three-leg forms and Lagrangians 



c) Three systems with the biquadratics pattern as on Figure 4c with two pairs of face 
equations being trapezoidal Hi, the third pair of face equations being rhombic Hi, 
and both tetrahedron equations Hi (for k = 1,2,3). 

All these six-tuples can be embedded into the lattice Z 3 by reflecting the cubes, see 
details of this construction in |Bolll| . The Lagrangian formulation of these systems 
continued to the whole lattice constitutes the main subject of the present paper. 

4. Three-leg forms and Lagrangians 



As already demonstrated in [BSlO], the Lagrangian structures of quad-equations are 
closely related to their three-leg forms, which, in turn, can be derived from the bi- 
quadratics related to the quad-equations. 

Lemma 4.1. For any of the polynomials Q k , Hi {k = 1,2,3) there exits a re-scaling of 
the function Q such that the re-scaled biquadratics depend only on one parameter, namely 
the parameter assigned to the corresponding edge: 

Q 1,2 (xi,x 2 ;ai,a 2 ) = a 2 )h 1,2 {xi, x 2 ; on), 

Q 1 ' 4 {xi,X4;ai,a 2 ) = ft(ai, a 2 )h 1,4 (xi, x 4 ; a 2 ), 
Q ' 3 (xi, x 3 ; a-i, a 2 ) = a 2 )/i 1,3 (xi, x 3 ; a\ - a 2 ). 

The proof consists of a direct verification with the scaling factors given in Sections 
A.l A. 2 In particular, for the polynomials Q k , due to their high symmetry, all three 



biquadratics are given essentially by one and the same non-degenerate symmetric poly- 
nomial h(x, y; a) = h(y, x; a): 

h 1 ' 2 (xi,x 2 ; a\) = h(xi,x 2 ; a{), /i 1,4 (xi, x 4 ; a 2 ) = —h(xi, x 4 ; a 2 ), 
/i 1,3 (xi, x 3 ; «i - a 2 ) = -h(xi,x 3 ; ct\ - a 2 ). 

For the polynomials H k , one still has 

h ' {x\, x 2 ; ax) = h(xi,x 2 ; «i), h 1 ' (x\, x^; a 2 ) = -h(x\, x 4 ; a 2 ), 

where the biquadratic polynomial h(x, y; a) is degenerate and non-symmetric. 

The signs on the right-hand sides of the latter equations are responsible for the choice 
of signs in the three-leg form of quad-equations: 

Lemma 4.2. Equations Q, and rhombic equations H k possess three-leg forms 
ip(xi,X2\ ai) - ip(xi,X4; a 2 ) - <p(xi,x 3 ; a\ - a 2 ) = 

and 

ijj(x 2 ,xi; ui) - V»(x 2 , x 3 ; a 2 ) - 4>{x 2 , x 4 ; ot\ - a 2 ) = 0, 
centered at x\ and at x 2 , respectively, with the functions vp,ip,(j) and (f) listed in Sec- 



tions A.l, A. 2 For equations Q e k , one has t/j = t/j = <f> = <j>. If x% = fi(Xi) denotes 
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4. Three-leg forms and Lagrangians 



the uniformizing changes of variables for yjr i {xi), where r l {xi) is the discriminant cor- 
responding to the vertex Xi, then the leg functions satisfy the following relations: 



ifj(x 1 ,x 2 ; ai) 
ip(x 1 ,x i ; a 2 ) 



dXi J h 1 < 2 (x 1 ,fr,a 1 ) 
f x * dU 



(x 1 ,x 3 ;a 1 - a 2 ) = wrp — / ttv i T> I 3 ) 

dXi J /i 1,d (xi,^ 3 ;ai - a 2 ) 



and 



ifj(x 2 ,xi; ai) 
rh(x 2 ,x 3 ; a 2 ) 



df 2 (x 2 ) ri d& 



ox 2 j /i 2 'i(x 2 ,a;ai; 

df 2 (X 2 ) [** d£ 3 



dX 2 J h 2 ' 3 (x 2 ,(, 3 ;a 2 ) 7 



(/){x 2 ,x 4 ;ai - a 2 ) = — — / -. (4) 



dX 2 J h 2 ' 4: (x 2 ,^; a\ - a 2 ) 
In addition, we have the following properties: 

• Symmetry: 

d d - 

d d 

-—-<f){xi y x 3 ;ai-a 2 ) = -——<j){x 3 ,xr,a.\-a 2 ), (6) 

d - d - 

-q^-J>{x 2 ,x 4 ; cti - ct 2 ) = g^4>{x4,x 2 ;a>i - a 2 ). (7) 

• Derivative with respect to the parameter: 



ip(x 1) X2;ai) = Q^rlogh ' (xi,X2;ai), (8) 

d - d 

Q^-ip(x 2 ,xi;ai) = ^^log/i 1,2 (xi,x 2 ;ai), (9) 

d d 

-4>(xi,x 3 ; ai - a 2 ) = -— — log h}' 3 (xi, x 3 ; a>i - a 2 ), (10) 



da\ ~ dX\ 

d - d 
-<f>(x 2 ,x 4 ; ct\ - ct 2 ) = --^rr\ogh 2A {x 2 ,x 4 ;ai - a 2 ). (11) 



da\ ' dX 2 

Proof. Similar statements for equations of the ABS list were already given in (ABS03, 
BS10J. Here we provide the reader with slightly more explanations why these statements 
hold true. Namely, the first statement follows from the result (due to V. Adler and 
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4. Three-leg forms and Lagrangians 



published as Exercise 6.16 in [BS08, p. 281]) that the three-leg form of an arbitrary 
multi-affine equation Q = 0, centered at x\, reads 



<3 1 ' 2 (xi,?2) J Q^(n,&) ; Q'- 3 (xi,U 



Due to Lemma 4.1 for equations H kl Q e k one can replace here Q l,:) by the biquadratics 
h 1 ' 3 depending on one parameter each. Furthermore, all three leg functions in the latter 
formula can be simultaneously multiplied by an arbitrary function of x\. The choice of 
this function as df\(X\)/dX\ ensures the symmetry properties of the derivative of the 
leg function with respect to the second argument: 

d dfxjXx) df 2 (X 2 ) 1 d - 

-^(xi,x 2 ;ai) = — — — — ■ -= „. - - , = — if; (x 2 , xr, ai) . 



dX 2 YK " " LJ 8X X dX 2 fcl. 2 (xi,x 2 ;ai) dXi 

The both diagonal biquadratics are symmetric but different, so that one has two different 
long leg functions (ft and (ft. The formulas for the derivatives of the leg functions with 
respect to parameters are checked by a direct computation. □ 

Remark. For trapezoidal equations H^, the notion of short and long legs might be some- 
what misleading. We will assume that, for both the rhombic and the trapezoidal ver- 
sions of Hi, the functions (ft and (ft correspond to two non-degenerate biquadratics (black 
and white, respectively), while the functions ift and ift correspond to four degenerate bi- 
quadratics, with the first argument (the center point of the three-leg form) being black, 
resp. white. Thus, the three-leg forms of a trapezoidal equation Hf,, as depicted on 



Figure lc read: 

4>(xi,x 2 ; ai) - ift(xi, x 4 ; a 2 ) + ift(xi,x^; a 2 - a\) = 
(centered at a black vertex x±), resp. 

(ft(x 3 , x 4 ; ax) - ift(x 3 , x 2 ; a 2 ) + ip(x 3 , x\\a 2 - a\) = 
(centered at a white vertex x 3 ). 

Remark. We recall that the three-leg forms of quad-equations yield the so called Laplace 
type equations of the stars of the vertices of the underlying quad-graph. These Laplace 
type equations say that the sum of long-leg functions for diagonals of all quadrilaterals 
adjacent to a vertex is equal to zero. For the rhombic equations the long- leg functions 
(ft, (ft are the same as for the type Q equations, so that they do not lead to new Laplace 
type equations. For the trapezoidal equations H^, the long- leg functions are ift, ift, so that 
one arrives at novel Laplace type equations. For instance, the Laplace type equations at 
the vertices X3 (black) and xq (white) on Figure [5] read: 

tp(x 3 , x 6 ; a 4 - a 3 ) - ift(x 3 , x 5 ; a 4 - «i) + ip(x 3 , xr, a 2 - aa) - 4>(x 3 , x 2 ; a 2 - 03) = 0, 
resp. 

■0(x 6 , x 8 ; a 6 - a 5 ) - tft(x 6 , xr, a 6 - a 3 ) + V5(x 6 , x 3 ; a 4 - a 3 ) - tft(x 6 , x 4 ; a 4 - a 5 ) = 0. 
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4. Three-leg forms and Lagrangians 




We are now in a position to construct functions which will play the role of Lagrangians 
in the action functionals for system of quad-equations. 

Lemma 4.3. For equations Q e k , H^, there exist functions L{X\,X 2 ;ai) (defined on 
directed edges starting at a black vertex and ending at a white vertex), 

A(Xi,X 3 ;ai - a 2 ) = A{X 3n X\]at - a 2 ) 

(defined on edges connecting two black vertices), and 

A (X 2 ,X 4i ; a\ - a 2 ) = A (X 4 , X 2 ; a.\ - a 2 ) 

(defined on edges connecting two white vertices) such that, upon changes of variables 
Xi = fi (Xi) introduced in Lemma there hold the following relations: 

d d 
tp (x 1 ,x 2 ;ai) = -T^L(X 1 ,X 2 ;ai) , tp (x 2 , x\\ a\) = q^-L (X 1 , X 2 ; ol{) , (12) 

d 

(j) (xi, x 3 ; ax - a 2 ) = ~t^ a {Xi,X 3 ; «i - a 2 ) , (13) 

d - 

(j)(x 2 ,x A ;ai -a 2 ) = -——A(X 2 ,X A ;a 1 - a 2 ) . (14) 

CA 2 

These functions (defined up to arbitrary additive functions of the respective parameters) 
can be chosen so that 

d 

L(Xi,X 2 ;ai) = log/i ' 2 (xi, x 2 ; ai) , (15) 



d i o 

A(Xi,X 3 ;ai - a 2 ) = log/i ' (xi, x 3 ; ai - a 2 ) , (16) 



da.\ 



Q 

A (A 2 , X4; a\ — a 2 ) = log h 2 ' 4 (x 2 , X4; qi — a 2 ) . (17) 



dai 
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4. Three-leg forms and Lagrangians 



The main theorem concerning the Lagrangians for a single quad-equation, which serves 
as the major ingredient of the proof of the flip invariance of the action functionals, reads: 

Theorem 4.4. For equations Hi, the Lagrangians L, A and A can be chosen so that the 
following relation holds on solutions of the equation Q(x\,x 2 , x 3 , X4) = 0: 



a) rhombic case with the enumeration of fields as on Figure lb 



L (X 1 ,X 2 ; at) + L (X 3 , X 4 ; a x ) - L {X U X A ; a 2 ) - L (X 3 , X 2 ; a 2 ) 

- A (Xi, X 3 ; ax - a 2 ) - A (X 2 , X 4 ; cti - ot 2 ) = 0; (18) 



b) trapezoidal case with the enumeration of fields as on Figure lc 



A(Xl, X 2 ; ai) + A(X 3 , X 4 ; a x ) - L{X X ,X±; a 2 ) - L(X 2 , X 3 ; a 2 ) 

+ L{X U X 3 ; a 2 - «i) + L(X 2 ,X 4 ; a 2 - ax) = 0. (19) 

For equations Q e k , the analogous statement holds true if one sets all Lagrangians L,A, 
and A equal (to A, say). 

Proof. The proof of this theorem is completely analogous to the one of Theorem 1 in 
[BS10] and is obtained by showing that the gradient of the function on the left-hand side 
vanishes on the solutions of Q = due to the equivalence of the equation Q = to its 
three-leg forms. □ 

It is important to notice that in a 3D consistent six-tuple of quad-equations possessing 
the tetrahedron property, the choice of Lagrangians of every quad-equation (be it the 
face equation or the tetrahedron equation) can be made in a consistent manner, i.e., on 
every edge (or diagonal) the two Lagrangians coming from the two different equations 
of the quadrilaterals sharing this edge coincide. This follows from the next result which 
was given for equations of the ABS list in |ABS09| . and was proven in full generality in 
|BoTTl] ): 

Lemma 4.5. Consider a 3D consistent six-tuple ([T]) of quad- equations possessing the 
tetrahedron property. Then: 

1. For any edge of the cube, the two biquadratics corresponding to this edge coincide 
up to a constant factor; for example 

£0,2 C 0,3 
= const, = const, -p^ = const. 

2. The product of these factors around one vertex is equal to —1; for example, for the 
vertex x one has: 

A ' 1 B ' 2 C°< 3 _ 
COA • ' PAZ ~ ~ L 

One can choose the normalization of the polynomials A, . . . ,C so that all such 
factors (for all 12 edges of the cube) are equal to —1. 
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5. Flip invariance of action functional 



As a consequence of this result, for consistent systems of quad-equations, each edge Z m 
carries a well defined leg function and a well defined Lagrangian, in the sense that these 
objects do not depend on the elementary square of Z m adjacent to this edge. Moreover, 
due to the construction of Lemma |3.2| the same is true for diagonals of elementary 
squares of Z m , which support leg functions and Lagrangians coming alternatively from 
edges of black and white tetrahedra supporting their own quad-equations. 

5. Flip invariance of action functional 

Following the main idea of |LN09| . we consider the functional for a multidimensionally 
consistent system of quad-equations defined on any quad-surfaces S in the Z m by the 
formula 

S= £ (20) 

where C((Tij) is the value of a discrete Lagrangian two-form on an elementary (oriented) 
square <7y = (n, n + e«, n + e% + ej, n + ej), 

£ (aij) = C (X,Xi,Xj;ai,aj) . 

For equations H k and Q e k defined on black-and-white (but not necessarily bipartite) 
quad-graphs, one can encounter eight types of elementary squares depicted on Figure [6] 

The corresponding values of the discrete two-form are: 



C(X, Xi,X 2 


Oil 


Oi2) 


= L(X,X 1 


Oil) 


-L(X,X 2 ; 


a 2 ) 


- A(Xi,X 2 ,a! 


- o 2 ); 


(a) 


C(X, Xi,X 2 


a\ 


Oi2) 


= L(X h X 


Oil) 


-L(X 2 ,X 


02) 


- A(Xi,X 2 ,ai 


- 02); 


(b) 


C(X, Xi,X2 




Oi2) 


= L(X, X x 


Oil) 


-A(X,X 2 ;o2) 


- L(X 2 ,Xi,ai 


- o 2 ); 


(c) 


C(X, X\,X2 


ai 


Oi2) 


= L(X U X 


Ot\) 


-A(X,X 2 


02) 


- L(Xi,X 2 ,ai 


- o 2 ); 


(d) 


C(X, Xi,X2 


Cti 


Oi2) 


= A(X,X 1 


Oil) 


-L(X,X 2 ; 


a 2 ) + L(Xi,X 2 ,a 2 


- 01); 


(e) 


C(X, X\,X2 


ai 


Oi2) 


= A(X,X 1 


Ol) 


-L(X 2 ,X 


02) 


+ L(X 2 ,Xi,a 2 


- 01); 


(f) 


C(X, Xi,X2 


Cti 


Oi2) 


= A(X,X 1 


Oi\) 


-A(X,X 2 ; 


02) 


-A(X u X 2 ,ai 


- o 2 ); 


(g) 


C(X, Xi,X2 


Oil 


Oi2) 


= A(X,X 1 


Ol) 


-A(X,X 2 ; 


02) 


-A(X u X 2 ,ai 


- 02). 


GO 



The choice of the sign of the third term in the expressions (c)-(f) is dictated by the 
requirement that the two terms with the Lagrangian L have opposite signs. 

With these definitions, one has the value of the action functional on any quad-surface 
in Z m , including Z 2 . In the latter case the action can be represented as the sum of La- 
grangians over all edges of Z 2 and over all diagonals running from north-west to south-east 
(or in the opposite direction, since they could be directed). The term "action functional" 
is justified by the following statement. 

Theorem 5.1. For any equation (rhombic or trapezoidal) or Q e k on the regular square 
lattice Z 2 7 the solutions are critical points of the action functional. 
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5. Flip invariance of action functional 




5. Flip invariance of action functional 



Our main result reads as follows: 
Theorem 5.2. For any 3D consistent system consisting of quad- equations H^. and/or 



Q e k (listed in Theorem 3.3), the action functional is flip-invariant in the sense that the 



following relation is satisfied on solutions: 

Ai£ {X, X 2 , X 3 ; a 2 , a 3 ) + A 2 C (X, X 3 ,Xx;a 3 , ax) + A 3 C (X, X 1 ,X 2 ; a u a 2 ) = 0. (21) 
Here Aj denotes the difference operator that acts on vertex functions according to 

AJ (x) = f (x t ) - f (x) . 
Proof We give the proof which generalizes the first proof of the corresponding result for 



equations from the ABS list given in [BSlOj. One adds the formulae from Theorem 4.4 



for the three faces and one tetrahedron adjacent to the point x. We give details for the 
three elementary cubes illustrated on Figure [4] and omit the details for the cubes obtained 
by the color inversion and rotations. 
Six-tuple from Figure 



Add the four relations 



L (X, X x ;a x ) + L (X 12 , X 2 ; ax) - L (X, X 2 ; a 2 ) - L (X 12 , X r , a 2 ) 

- A (X, Xi 2 ; cl\ - a 2 ) - A (Xx,X 2 ; oi\ - a 2 ) = 0, 

- L (X, X r , ai ) - L (X 13 ,X 3 ; Ql ) + L (X, X 3 ; a 3 ) + L (A 13 , X x - a 3 ) 

+ A (X, X 13 ; ai - a 3 ) + A (X x , X 3 ; ax - a 3 ) = 0, 

L (X, X 2 ; a 2 ) + L {X 23 , X 3 ; a 2 ) - L (X, X 3 ; a 3 ) - L (X 23 , X 2 ; a 3 ) 

- A (X, X 23 ; a 2 - a 3 ) - A (X 2 , X 3 ; a 2 - a 3 ) = 0, 



A (X, X 23 ; a 2 - a 3 )+A (X 12 , X 13 ; a 2 - a 3 )-A (X, X 13 ; ax - a 3 )-A (X l2 , X 23 ; ax - a 3 ) 

+ A (X, X l2 ; ax - a 2 ) + A (X 13 , X 23 ; ax - a 2 ) = 0. 

The result reads: 

L(X 12 ,Xx;a 2 ) + L{X 23 ,X 2 ;a 3 ) + L(X 13 ,X 3 ;ax) 

-L (Xx 3 , Xx; a 3 ) - L (X 12 , X 2 ; ax) - L (X 23 ,X 3 ; a 2 ) 

-A (Xx2,Xx 3 ; a 2 - a 3 ) - A (Xx 2 ,X 23 \ a 3 - ax) - A (X l3 , X 23 ; ax - a 2 ) 

+A(X 2 ,X 3 ;a 2 - a 3 ) + A (Xx,X 3 ; a 3 — ax) + A (Xx,X 2 ; ax - a 2 ) = 0. 



which is the flip invariance relation (21 ) in this case. 



Six-tuple from Figure J^b: Add the four relations 



A (X, Xx; ax) + A (X 2 , X x % ax) - A (X, X 2 ; a 2 ) - A (X X ,X 12 ; a 2 ) 

— A (X, Xx2', ax — a 2 ) — A (Xx,X 2 ; ax — a 2 ) = 0, 
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5. Flip invariance of action functional 



- A (X, Xl; ax) - A (A" 3 , A 13 ; a ± ) + L (X , X 3 ; a 3 ) + L (X x , X l3 ; a 3 ) 

+ L (A, A i3 ;ai - a 3 ) + L (Xi,X 3 ;a!i - a 3 ) = 0, 

A (A, A 2 ; a 2 ) + A (A 3 , A 23 ; a 2 ) - L (A, A 3 ; a 3 ) - L (A 2 , A 23 ; a 3 ) 

- L (A, A 23 ; a 2 — a 3 ) — L (A 2 , A 3 ; a 2 - a 3 ) = 0, 

L (A, A 23 ; a 2 - a 3 )+L (A i2 , A i3 ; a 2 - a 3 )-L (A, A i3 ; ai - a 3 )-L (A i2 , A 23 ; ai - a 3 ) 

+ A (A, A i2 ; ai - a 2 ) + A (A i3 , A 23 ; a x - a 2 ) = 0. 

The result reads: 

A(Ai, Ai 2 ;a 2 ) + L(A 2 , A 23 ;a 3 ) + A(A 3 ,Ai 3 ;ai) 

-L(Xi, Ai 3 ;a 3 ) - A (A 2 , A i2 ; ai) - A (A 3 , A 23 ; a 2 ) 

-L (X12, X 13 ; a 2 - a 3 ) + L (Ai 2 , A 23 ; ol\ - q 3 ) - A (A i3 , A 23 ; a\ - a 2 ) 

+L(A 2 , A 3 ;a 2 - a 3 ) - L(X 1 ,X 3 ;ai - a 3 ) +A(Ai,A 2 ;ai - a 2 ) = 0. 



which is the flip invariance relation (21 ) in this case. 



Six-tuple from Figure 4c: Add the four relations 



L (A, X r ,a x ) + L (A 12 , A 2 ; c*i) - L (A, A 2 ; a 2 ) - L (A 12 , A l5 a 2 ) 

- A (A, Ai 2 ;ai - a 2 ) - A(Ai, A 2 ;ai - a 2 ) = 0, 

- L (A, A i; ai ) - L (A 3 , A 13 ; ax) + A (A, A 3 ; a 3 ) + A (A 1; A 13 ; a 3 ) 

+ L (A, A i3 ;ai - a 3 ) + L (X 3 , Ai; ot\ - a 3 ) = 0, 

L (A, A 2 ; a 2 ) + L (A 3 , A 23 ; a 2 ) - A (A, A 3 ; a 3 ) - A (A 2 , A 23 ; a 3 ) 

- L (A, A 23 ; a 2 - a 3 ) - L (A 3 , A 2 ; a 2 - a 3 ) = 0, 

L (A, A 23 ; a 2 - a 3 )+L (A i2 , A i3 ; a 2 - a 3 )-L (A, Ai 3 ; ai - a 3 )-L (Ai 2 , A 23 ; a>i - a 3 ) 

+ A (A, A i2 ; ai - a 2 ) + A (A i3 , A 23 ; a\ - a 2 ) = 0. 

The result reads: 

L(Ai 2 , Ai;a 2 ) + A(A 2 ,A 23 ;a 3 ) + L(A 3 , Ai 3 ;ai) 

-A (Ai, Ai 3 ; a 3 ) - L (A i2 , A 2 ; a x ) - L (A 3 , A 23 ; a 2 ) 

-L (A12, A i3 ; a 2 - a 3 ) + L (Ai 2 , A 23 ; ai - a 3 ) - A (A i3 , A 23 ; ai - a 2 ) 

+L(A 3 , A 2 ;a 2 - a 3 ) - L(A 3 , X x \a x - a 3 ) + A(Ai, A 2 ;ai - a 2 ) = 0. 



which is the flip invariance relation (21 ) in this case. 



□ 
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A. List of Quad-Equations 
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A. List of Quad-Equations 
A.l. Type Q 

Q{ : Q = cci (x\x 2 + X3X4) - a 2 (xix 4 + x 2 x 3 ) - (cci - a 2 ) (xi%3 + x 2 x 4 ) 
+e 2 aiQ 2 (ai - a 2 ), 
X{ — X{ , 

K(ai,a 2 ) = 2aia 2 (ai - a 2 ), 

1 / {x\ - x 2 f 



^2;aiJ = 7: I : e 2 ai 

2 V ai 

( 1 . A" 4 - A" 2 + eai 



(/>(xi,x 2 ;ai) = < 



log 7 7 ' — , ife/O, 
e Ai — A 2 — ecci 

v v ' ife = 0; 

I, Ai — A 2 



Q|: Q = ai {x\x 2 + X3X4) - a 2 (xix 4 + x 2 x 3 ) - (ai - a 2 ) (xix 3 + x 2 x 4 ) 
+aia 2 (ai - a 2 ) + eaia 2 (ai - a 2 ) (xi + x 2 + x 3 + x 4 ) 
-e 2 «ia 2 (ai - a 2 ) (a 2 - aia 2 + a^), 

/X(l + eX), ife^O, 
Xi = < 

]X, if e = 0, 

a 2 ) = 2aia 2 («i — a 2 ), 
/i 1,2 (xi, x 2 ; ai) = ((xi - x 2 ) 2 - af - 2eaf (xi + x 2 ) + e 2 af ) , 
(Xi - X 2 + ai) (1 + e(X 1 + X 2 + ai)) 



(p(x 1 ,x 2 ;a 1 ) 



log 
log 



(Xi - X 2 - ai) (1 + e(Xi + X 2 - ai)) ' 
Ai - X 2 + ai 



A~i — X2 — ai ' 

Q\: Q = sinh(2ai) (xix 2 + x 3 x 4 ) — sinh(2a 2 ) (x\x 4 + x 2 x 3 ) 
- sinh(2(ai - a 2 )) (xix 3 + x 2 x 4 ) 
—A5 2 e 2 sinh(2ai) sinh(2a 2 ) sinh(2(ai - a 2 )), 

f2(5ecosh(2X), if e / 0, 
Xl ~\exp(2X), if e = 0, 

n(a±, a 2 ) = sinh(2ai) sinh(2a 2 ) sinh(2(ai — a 2 )), 



if e / 0, 
if e = 0; 
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A. List of Quad-Equations 



12; ( e^xi - e~ ai x 2 ) (e-^ai - e 2a ^x 2 ) | ( , 22 

sinh(2ai) 




+ 45 2 e 2 sinh(2ai), 



sinh(Xi + X 2 + ai) sinh (X x - X 2 + «i) /n 
lQ g ....... , v , v , • , . ^ ; \ ' lf e/0, 



log- 



sinh (X\ + X 2 — ai) sinh (Xi — X 2 — a{) 
sinh (Xi — X 2 + «i) 



sinh (Xi — X 2 — ai) 



+ 2qi, 



if e = 0. 



A.2. Type H 4 

if^: Q = (xi — x 3 ) (x 2 - x 4 ) + (a 2 - ai) (l - e 2 x 2 X4). 

— X^, 

K(ai,a 2 ) = 2(ai — 02), 

1 - e 2 x 2 

h 1 ' 2 (xi,x 2 ;ai) = — -, 



h 1 ' 3 (xi,x 3 ;ai - a 2 ) 



h 2 ' A (x 2 , X4; ai — a 2 
ip (xi,x 2 ;ai) = 



1 / (xi - x 3 ) z 

2 I ai — «2 

1 / (x 2 - x 4 ) z 

2 \ ai — a 2 



e 2 (ai - a 2 ) 



-log- =r, lfe^O, 

e 1 - eX 2 

2X 2 , if e / 0, 

Xi - e 2 aiX 2 



ip (x 2 , xi; ai) 



1 - e 2 X 2 



2Xi 



-, ife/0, 
if e = 0, 



</>(xi,x 3 ;ai - 02) 
4> (x 2 , X4; ai - 02) = 2 



f 1 Xx - X 3 + e (ai - a 2 ) if £ _^ Q 
e Xi — X 3 — e (ai — a 2 ) ' 



cti — a 2 
, Xi — X 3 
ai — a 2 
X 2 — X4 ' 



if e = 0, 



H%: Q = (xi — x 3 ) (x 2 - x 4 ) + (a 2 - «i) (xi + x 2 + x 3 + x 4 ) + a\ — a\ 
-| («2 - «i) (2x 2 + "i - 
(X(l + eX), ife/0 



(a 2 - ai) (2x 2 + «i + 0S2) (2x 4 + «i + a 2 ) - - (a 2 - a>i) , 



[Xi if e = 

k(o>i,o> 2 ) = 2(ai — a 2 ), 

h 1 ' 2 (xi, x 2 ; ai) = xi + X2 + cti - e (x 2 + ai) 2 



(i = 1,3), x, = X (i = 2,4), 



/i 1,3 (xi,x 3 ;ai - a 2 ) 



2 (ai - a 2 ) 
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A. List of Quad-Equations 



■ ((xi - X3) 2 - (ai - a 2 ) 2 - 2e (ai - a 2 ) 2 (zi + £3) + e 2 (ai - a 2 ) 4 ) , 



/T' 4 (x 2 , x 4 ; ai - a 2 ) = - 



2 (ai - a 2 ) 
Xi + X 2 + ai 



■0 (xi,x 2 ;ai) : 
ijj (x 2 ,xi;ai) : 

(j)(xi,x 3 ;a>i - a 2 ) 



(x 2 - x 4 ) 2 - (ai - a 2 ) 2 ). 
, if e / 0, 



1 ° g l + e(X 1 -X 2 -a 1 ) 
log (Xi + X 2 + ai) , if e = 0, 

log((Xi + X 2 + ai) (1 + e(Xi - X 2 - ai))), if e / 0, 
log(Xi + X 2 + ai), if e = 0, 

(Xi - X 3 + ai - a 2 ) (1 + e(Xi + X 3 + ai - a 2 )) ., , . 
log^ Az ef, if e 7^0, 



log 



(Xi - X 3 - ai + a 2 ) (1 + e(Xi + X 3 - a x + a 2 )) ' 



X3 + ai — a 2 



0(x 2 ,x 4 ;ai - a 2 ) = log 



X\ — X 3 — ol\ + a 2 
X 2 — X4 + a\ — a 2 



if e = 0, 



e 2 x 2 X4 \ 



X 2 — X/l — ai + a 2 ' 
Q = e 2ai (xix 2 + x 3 x 4 ) - e 2 " 2 (x : x 4 + x 2 x 3 ) + (e 4ai - e 4a?2 ) f <5 2 + - ; 
f2<5ecosh(2Xi), if e / 



\exp(2Xi), if e = 

K(ai,a 2 ) = (e 4Ql - e 4a2 ) /2, 

/i 1 ' 2 (xi, x 2 ; ai) = -2 ( xix 2 + 5 2 e 2a?1 + 6 '' 



(i = 1,3), x l = exp(2X i ), (i = 2,4), 



2^2 
2 

2ai 



/i 1,3 (xi,x 3 ;ai - a 2 ) 



( e ai-«2 Xl _ e "2-ai X3 ) ( e «2-ai Xl _ e «i-a2 X3 ) 



sinh(2(ai — a 2 )) 



- A5 2 e 2 sinh(2(ai - a 2 )), 



h 2 ' A (x 2 ,x 4 ;ai - a 2 ) 



(e ai - a2 x 2 - e a2 ~ ai x 4 ) (e a2 " ai x 2 - e ai " a2 x 4 ) 



ip (xi,x 2 ;ai) = < 



ip (x 2 , xi; ai) = < 



log 

log 
log 

log 



2X2-2X1 1 J e 2ai 



ee 2Xi+2X 2 _|_ (5 e 2ai 
g2A"i+2X2 _|_ ^2g2cii 



sinh(2(ai — 0:2)) 
, ife/0, 

, if e = 0, 



(ee 2X 2-2X 1 + §e 2a^ ( ee 2X 1+ 2X 2 + de 2 ai ) 
p 2ai 



e 2Xi+2X 2 _|_ ,52 e 2oi ' 



, if € ^0, 
if e = 0, 
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(f>(xi,x 3 ;a>i - a 2 ) = < 



^ sinh (Xi + X 3 + ai - a 2 ) sinh (Xi - A 3 + ct\ - a 2 ) ^ e _^ q 
sinh (X\ + X3 — ai + a 2 ) sinh (X\ — X3 — ct\ + a 2 ) ' 

sinh(Ai - X 3 + ai - a 2 ) 

log ■ , / y : r + 2(ai-a 2 ), if e = 0, 

smn(Ai — A3 — at\ + a 2 ) 



- sinh(A' 2 - X 4 + ai - a 2 ) 

(x 2 , x 4 ; ai - a 2 ) = log . ■ + 2(a x - a 2 ). 

smh(A 2 — A4 — ai + a 2 ) 

Remark. The transition from Q 3 to Q3 and from to ffg is performed by first shifting 
Xi — > Xi — log((5e) (for i = 1, 2, 3, 4 for equation Q| and for i = 1, 3 for equation 
and then sending e — >■ 0. 
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